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Abstract: The prospects to test the hypothesis of intrinsic charm quarks in the 
, proton are investigated. We consider how this component can be directly or indi- 

J> \ rectly probed in deep inelastic scattering at HERA and in fixed target experiments 

and find that an overlooked signal might be present in existing NMC data. Applying 
. the intrinsic charm model to hadron collisions we compare the resulting charm pro- 

duction cross-sections with those based on standard perturbative QCD and available 
data. Extrapolating to higher energies we obtain predictions for charm production 
■ at the Tevatron and LHC. 

*5 : 

^ 1 Introduction 

The hypothesis of intrinsic charm quarks in the proton was introduced []l]] as an attempt to 
understand a large discrepancy between early charm hadroproduction data and leading order 
perturbative QCD (pQCD) calculations. The data from ISR (see refs. in Q) were about an 
order of magnitude higher than the prediction and had a rather flat distribution in longitudinal 
momentum compared to the sharp decrease with Feynman xp = P\\/p m ax expected from pQCD. 
As more data have been collected at different energies and next-to-leading order (NLO) pQCD 
calculations have been made, the discrepancy has largely disappeared. Still, however, there are 
certain aspects of the data which are difficult to understand within the pQCD framework, but 
are natural if the intrinsic charm hypothesis is basically correct. This concerns the mentioned xf 
distributions of leading charm in hadroproduction |||, the dependence on the nuclear number 
(A a ) of J/tp production Q and double J/ift production ||. 

Thus, although the main features of charm production can now be understood by conven- 
tional pQCD, certain aspects of the data indicate the presence of some additional mechanism. 
Intrinsic charm may be such a process, which gives a small contribution to the inclusive cross 
section but could dominate in some regions of phase space. 

The hypothesis of intrinsic charm (IC) amounts to assuming the existence of a cc-pair as 
a non-perturbative component in the bound state nucleon [pj. This means that the Fock-state 
decomposition of the proton wave function, \p) = a\uud) + (3\uudcc) + contains a small, 
but finite, probability /3 2 for such an intrinsic quark-antiquark pair. This should be viewed as a 
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quantum fluctuation of the proton state. The normalization of the heavy quark Fock component 
is the key unknown, although it should decrease as I/tuq. Originally, a 1% probability for 
intrinsic charm was assumed, but later investigations indicate a somewhat smaller but non- 
vanishing level; ~ 0.3% f§ and (0.86 ± 0.60)% §. 

Viewed in an infinite momentum frame, all non-perturbative and thereby long-lived compo- 
nents must move with essentially the same velocity in order that the proton can 'stay together'. 
The larger mass of the charmed quarks then implies that they take a larger fraction of the proton 
momentum. This can be quantified by applying old-fashioned perturbation theory to obtain the 
momentum distribution [ffll 



P(p — » uudcc) oc 



in 



E 



m 



2 



(1) 



in terms of the fractional momenta Xi of the five partons i in the \uudcc) state. The probability 
P 2 of this state is related to the normalisation N$ of ref. Neglecting the transverse masses 
{m\ = m? + Pj_) of the light quarks in comparison to the charm quark mass results in the 
momentum distribution 

P(x 1 ,X 2 ,X 3 ,X c ,Xc) OC ; C 5(1 - X! - X 2 ~ X 3 - X c - X Z ) (2) 

where the (transverse) charm quark mass is absorbed into the overall normalisation. This 
function favours large charm quark momenta as anticipated. In fact, one obtains (x c ) = 2/7 by 
integrating out the remaining degrees of freedom x,. 

An intrinsic cc quantum fluctuation can be realised through an interaction, such that charmed 
particles are created. In proton-proton collisions this could certainly happen through a hard 
interaction with such a charm quark, but the cross section is then suppressed both by the small 
probability of the fluctuation itself and by the smallness of the perturbative QCD cross section. 
The charm quarks may, however, also be put on shell through non-perturbative interactions that 
are not as strongly suppressed ||]. Deep inelastic lepton-proton scattering would be another 
possibility, where this intrinsic charm component could be more directly probed. To investigate 
these possibilities we have constructed a model based on refs. [[[], ^ ||]. In part we use Monte 
Carlo techniques to simulate explicit events giving a powerful method to extract information on 
various differential distributions. 

One may also consider the extension from intrinsic charm to intrinsic bottom quarks. Since 
the overall probability scales with 1/mQ, one expects all cross sections for intrinsic bottom to 
be about a factor ten lower. The x-distributions for bottom quarks should be somewhat harder 
than for charm, but since light quark masses have already been neglected in Eq. (||) it also applies 
for intrinsic bottom and the same distribution as for charm is obtained in this approximation. 

In section 2 we investigate IC in deep inelastic scattering (DIS) in fixed target and collider 
mode. In the former, we explicitly consider the experiment New Muon Collaboration (NMC) 
at CERN and find that an interesting number of intrinsic charm events could be present in 



their data samples. For ep collisions in HERA, we extend a previous study [ 10 1 to include 
the possibility of scattering on a light quark such that the intrinsic charm cc-pair is freed and 
give charmed particles in the forward proton remnant direction. The case of hadronic charm 
production is studied in section 3. The measured charm production cross sections are compared 
to those calculated in pQCD to constrain the allowed amount of IC. Here, an important issue 
is the energy dependence of the cross section for IC, where two alternatives are considered. 
Based on this, we provide differential distributions of charm production from IC in comparison 
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Figure 1: Deep inelastic scattering on a proton containing an intrinsic cc-pair; a) scattering on 
the charm quark, b) scattering on a light valence quark. 



to the standard pQCD treatment, as applied to the Tevatron and LHC. Finally, section 4 gives 
a concluding discussion. 

2 Intrinsic charm in deep inelastic scattering 
2.1 General framework 

In DIS it should be possible to measure the effective IC parton density by direct scattering on an 
intrinsic (anti)charm quark (Fig.[l]a). The parton density is obtained from Eq. @ by integrating 
out all degrees of freedom except the momentum fraction x of the charm quark (or antiquark), 
resulting in 

c IC (x) = /3 2 1800x 2 j^(l-x)(l + 10x + x 2 ) + 2x(l + x)lnxj (3) 

The DIS charm cross section is given by 
d 2 <7 2ira r 



dxdQ 2 xQ 4 



l + il-yf FZ(x,Q 2 ) (4) 



in terms of the conventional DIS variables; Bjorken x = Q 2 /2P • q, y = P ■ qj P ■ pi where 
P,Pe is the momentum of the initial proton and lepton, and q is the four momentum transfer of 
the exchanged photon. Only photon exchange contribution is here included since the high-Q 2 
region with Z° exchange contributions can be neglected in our applications. The general charm 
structure function is 

F 2 c (x, Q 2 ) = e 2 [xc(x, Q 2 ) + xc(x, Q 2 )} (5) 

and F\ is obtained if the intrinsic charm quark density cjc{x, Q 2 ) is used. The Q 2 dependence 
from normal leading log GLAP equations have been calculated for IC in ||], but can be taken 
into account through a simple extension of the parameterisation in Eq. (||) ]T3] . 

We have included an option to treat DIS on intrinsic heavy quarks in the Monte Carlo 



program Lepto [11] using this formalism and are therefore able to simulate complete events 
based on this IC model. As for normal DIS events, higher order perturbative parton emissions 
from the incoming and scattered quark are included through parton showers. Hadronisation 



and particle decays are then performed according to the Lund model [12] in its Monte Carlo 



implementation Jetset |L3| 
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In addition to this direct scattering on an intrinsic charm quark, one may also consider the 
case of DIS on a light (valence) quark (Fig. |]b) such that the intrinsic cc-pair in the proton is 
'liberated' and give charmed particles in the phase space region corresponding to the 'spectator' 
proton remnant. To simulate this case also, we have made add-on routines to Lepto based on 
the following simple model. The normal electroweak DIS cross section is first used to treat the 
basic scattering on a light quark. The corresponding quark density function is here obtained 
from Eq. (^) by integrating out all degrees of freedom, except x\ for the light (valence) quark to 
scatter on, 

qic(xi) =/? 2 6(l-xi) 5 . (6) 

The probability /3 2 of the \uudcc) state is included here and gives the normalisation of the cross 
section for this intrinsic charm process. 

From the remaining qqcc system, one then considers the production of Z)-mesons, A c and 
J/tp through a coalescence [Q, [| of the relevant partons in Fig.|l|b. The longitudinal momentum 
distributions for these charmed particles are obtained, for a given x\, by integrating over the 
other momentum fractions in Eq. (g) and inserting a <5-function for the particle to be produced, 
i.e. 5{xQ — x q — x s ), 6(xa c — x u — Xd — x c ) or 5(xj/^ — x c — x 5 ). Considering momentum 
fractions relative to the proton remnant, i.e. after x\ is removed, corresponds to the variable 
substitution x — » £ = x/(l — x\). This gives expressions that factorise into an aq-dependent part 
and a ^-dependent part. Since x\ is already chosen, the former part can be absorbed into the 
normalisation of the ^-distribution. Thus, normalising to overall unit probability, one obtains 
the probability distributions 
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P B 300 j£ 4 hi£ + (1-O 4 ln(l-0 + 
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P Ac 300 (l-O 2 {e[6 -5£] +2[i 2 + 4£ + 3] ln(l - £)} (8) 
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dN W _ n on , , < ,< -■', 

an 



where Pjj,P\ c and Pj/^p are the relative probabilities to form the three different charmed par- 
ticles. These distributions are illustrated in Fig.|2| 

The overall normalisation is already accounted for (by /3 2 ) in the DIS cross section and only 
the relative probabilities Pq, Pa c an d Pj/^> need to be specified. Following the investigations 
in |Q, [| for the similar case in hadroproduction, we take 66% D, 33% A c /S c and 1% J/ifj. In 
each case there is a remaining parton in the proton remnant, i.e. a qc diquark, a c quark and 
a qq valence diquark, respectively. These are all anti-triplets in colour and are connected with 
a Lund string to the scattered quark to form a singlet system that hadronises using the normal 
Lund model. Note, that the charm (anti)quark at the proton remnant end of the string also 
gives rise to a charmed hadron. Before hadronisation, however, the scattered quark may radiate 
partons as treated by parton showers in Lepto. 



In the following, we extend the previous study |1C] of direct DIS on IC in two ways. First, 
a more dedicated study of direct IC scattering in connection with the New Muon Collaboration 
(NMC) experiment is made. Secondly, we consider the above model for indirect IC scattering 
applied to ep collisions at HERA and compare with direct IC scattering as well as conventional 
boson-gluon fusion into cc from pQCD. 
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Figure 2: Distribution of rescaled longitudinal momentum £ = x/(l — x\) of (a) D-mesons, 
(b) A c and (c) J ftp produced through coalescence in the remnant system from a DIS scattering 
on a valence quark in the intrinsic charm \uudcc) state. The distributions shown in the scaled 
variable are normalised to unit integral as discussed in the text. 



2.2 Fixed target muon scattering (NMC) 

In case of fixed target experiments the interest is focused on the direct scattering on intrinsic 
charm. This gives a high energy charmed particle in the current direction. The resulting high 
energy muon from a semi-leptonic decay can then be used as a signal for charm. In the indirect 
IC process, however, the charmed particles will emerge in the target region and resulting muons 
will have low energy and hence be more difficult to identify. Therefore we will not consider this 
indirect mechanism in the fixed target case. 

We thus concentrate on the direct scattering on intrinsic charm quarks and use our Monte 
Carlo model to simulate events corresponding to the NMC experimental situation. This means 
having a 280 GeV muon beam on a stationary proton target and applying the cuts Q 2 > 2 GeV 2 
and W 2 > 100 GeV 2 . 

From the simulated intrinsic charm events we extract the structure function F% ', using 
Eq. (Q), resulting in the three curves in Fig. ||] for different intrinsic charm probabilities. For 
comparison, we have calculated F£ arising from the conventional boson-gluon fusion (BGF) 
process 7*5 — > cc. This was obtained from a similar Monte Carlo simulation using the Aroma 
program @, which is an implementation of this process using the leading order pQCD matrix 
elements with the proper charm quark mass effects included (together with parton showers and 
hadronisation as in Lepto). The large- x feature of IC is clearly seen and results in a dominance 
of IC over pQCD at large enough x. The cross-over point depends on the unknown absolute 
normalisation of IC as illustrated. 

One should here consider the experimental information on the inclusive F% available from the 
European Muon Collaboration (EMC). Their original analysis gave an upper limit for IC (/? 2 ) 



of 0.6% at 95% CL [15]. In a later analysis ||, the QCD evolution of the IC structure function, 
as mentioned above, and charm mass effects were taken into account giving positive evidence for 
IC at the level 0.3%. A very recent investigation [Q] (in parallel with our study) has improved 
the theoretical treatment by including next-to-leading order corrections both for the IC and 
pQCD processes. Based on the same EMC data they find that an intrinsic charm contribution 
of (0.86 ± 0.60)% is indicated (in the bin of large energy transfer with mean v = 168 GeV). 

Thus, there is some evidence for IC from the inclusive F£ measurement of EMC, but it is 
not conclusive. One therefore needs to consider whether further information can be obtained. 

The success of the experimental method to tag charm through muons depends very much on 
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Figure 3: The charm structure function F^x) obtained from our Monte Carlo simulation of DIS 
on intrinsic charm quarks, with the different indicated normalisations (f5 2 ), and the conventional 
pQCD boson gluon fusion process j*g — > cc (solid curve). The error bars indicate the potential 
statistical precision of the discussed NMC data sample. 



the amount of muons from other sources. The background to IC from perturbatively produced 
charm, as given above, is in a sense an irreducible background. There is, however, also a 
severe background from the decays of light hadrons, mostly pions and kaons, that are copiously 
produced. Since these particles have a much longer life-time than charm, their effect can be 
reduced by either having the muon detectors close to the target or having hadron absorbers to 
filter out the prompt muons. 

The normal experimental arrangement of NMC is to have little material between the tar- 
get and the detectors. This implies, however, a substantial muon background from tt and K 
decays. However, in some runs of the NMC experiment a set-up more suitable for a charm 
search was used. In the heavy target configurations during 1987-88 there were hadron absorbers 
immediately downstream of the targets. We have developed an analytic method to calculate the 
suppression of this muon background due to the hadron absorption. The details are reported in 
lilj , where it is also applied to this NMC set-up. Thus, we calculate the development of the flux 
of produced mesons within the targets, in the calorimeters/absorbers and in the empty space 
between these elements. The resulting muon flux from hadron decays was found to be reduced 
to only ~ 6% of the flux without such absorption. In addition, we found that the calorimeter 
between the target and the absorber could also be used as a target. Charmed hadrons can also 
be absorbed, but due to their much shorter life-time the effect is negligible. For example, the 
decay length of a is at most ~ 1/10 of its interaction length (corresponding to a momentum 
of 180 GeV, which is about the highest to be expected). 

An experimentally realistic study of signal and background muons has been performed using 
Monte Carlo simulations. The results are displayed in Fig. ||. The background from light hadron 
decays was obtained by applying the above suppression factor on the yield from simulated normal 
DIS events. Here, the standard version of Lepto was used in the same kinematic region as 
specified above. The muons from IC and the pQCD process are obtained with our modified 
Lepto version and Aroma, respectively. 

The resulting energy spectra of muons are shown in Fig.|I| for three sets of cuts. Fig.|]a is 
for the above mentioned standard cuts; Q 2 > 2 GeV to ensure DIS and W 2 > 100 GeV 2 to 
have enough energy in the hadronic system for charm pair production. Fig.||b have extra cuts 
motivated by the NMC experiment; the scattered muon in the energy range 42 < E^i < 224GeV 
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Figure 4: Differential cross section versus energy of muons + /x + ) from the decay of charmed 
particles originating from a 1% intrinsic charm component in the proton (solid curves), from 
pQCD 7*5 — ► cc (dashed curve) and from decay of light mesons (mainly n,K) (dotted curves); 
all in 280 GeV muon-proton scattering with cuts: (a) Q 2 > 2GeV and W 2 > 100 GeV 2 , (b) 
NMC cuts on the scattered muon 42 < E^i < 224 GeV and 9 > 10 mrad and on secondary 
muons Eu2 > lOGeV, (c) Eu\ > 100 GeV and xbj > 0.2 to enhance IC. An overall reduction 
for muons from tt and K decays to 6% has been applied to account for the discussed hadron 
absorption effect in NMC. 



and at angle 6 > 10 mrad to emerge from the beam and secondary muons (from other sources) 
having a large enough energy E^ > WGeV to be indemnified. In Fig.|4|c the additional cuts 
Epi > 100 GeV and xsj > 0.2 were also made to identify the highest energy muon as the 
scattered one and select the large- x region. This last requirement is ment to enhance intrinsic 
charm events over the background processes (cf. Fig. ||). As demonstrated, it is then possible 
get muons from IC to the same level as the backgrounds. 

With the additional cuts the absolute rate is lowered such that large statistics data samples 
are needed. This may be satisfied by the NMC run in 1988 when data were collected for 
9.5 • 10 11 muons of 280 GeV on four targets of thickness 145 g/cm 2 each, which corresponds to 
an integrated luminousity of 334pb~ 1 . The statistical precision that can be expected in F% for 
this sample is illustrated in Fig. || for a few x-bins. This should give ~ 200 events with muons 

+ of energy larger than 10 Gel/ and fulfilling the stricter set of cuts, i.e. corresponding 
to the integrated distribution in Fig.|3|c. In the same way one obtains ~ 85 muons from charm 
produced by the pQCD process and ~ 340 muons from tt, K decays after the hadron absorpion 
reduction as estimated with the above mentioned analytic method. Additional data based on 
8.4 • 10 11 muons of 200 GeV energy could also be analysed in this way. 

In order to reduce the background from 7r, K decays further, one might also use secondary 
vertices based on the difference in decay length of light mesons as compared to charmed particles. 
We have investigated [16] this by considering the impact parameter of the muon tracks with 
respect to the primary interaction vertex. The impact parameter for muons from charm is 
always less than 0.5 cm, whereas it is typically larger for those from light mesons. However, the 
hadron absorbers that reduce the rate of muons from light mesons, also give rise to multiple 
scattering of the muons passing through. This causes a smearing of the reconstructed impact 
parameter which we find [l(| has a Gaussian distribution of width 3.5cm. Thus, the measureable 
impact parameter of muons from charm does not seem to provide a useable signature, although 
a safe conclusion would require a detailed simulation of the experimental conditions. 

In these estimates for IC we have assumed a 1% normalisation {(3 2 ) of the IC component 
and leave a trivial rescaling to any smaller value that may be preferred. Based on the above, 
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we conclude that it seems possible to find an intrinsic charm signal down to the 0.2% level in 
the NMC data through an excess of muons as compared to the expected background (e.g. an 
excess of about 40 over a background of 425 muons). The same sensitivity can be estimated 
from the precisions of F$ as shown in Fig.||. The determination of the BGF background is here 
essential, but uncertainties in the overall normalisation of the theoretical calculation can be well 
determined at small- x where this process dominates F% and the statistical precision is high. 

One should note that at large Q 2 , the intrinsic and extrinsic (BGF) charm processes have 
different scaling behaviour with the charm quark mass; 1/m 2 and log(m 2 ), respectively. This 
can give an additional handle to separate them in a data sample. 



2.3 Electron-proton collider 

In HERA at DESY, 30 GeV electrons (or positrons) are collided with 820 GeV protons. The di- 



rect DIS on intrinsic charm quarks in the proton has been investigated previously [10|. Since the 
intrinsic charm quark typically has a large x, it will be scattered at a rather small forward angle 
close to the proton beam direction. This makes it difficult to detect muons from semileptonic 



decays and only a fraction of the events are observable. For details we refer to [ 10 [ . 

The above discussed case of indirect intrinsic charm scattering may also occur. The charm 
particles are then produced in the forward-moving proton remnant system and will emerge in 
the very forward direction and therefore be hard to detect. Still, it is worthwile to estimate the 
rate and the distribution in pseudo-rapidity of the emerging muons. 

In a recent paper 17] these different IC processes in ep collisions were discussed at a theoret- 
ical and qualitative level. Here, we concentrate on a phenomenological study of these processes 




Figure 5: Distribution in pseudorapidity (n = —In tan(8/2) ) for muons from semileptonic decays 
of charmed particles produced through different DIS processes at HERA: scattering on a light 
(valence) quark in the intrinsic charm proton state \uudcc) (solid curves), direct scattering on 
an intrinsic charm state (dashed curves) and the conventional pQCD j*g — > cc process (dotted 
curves). In each case, the upper curve correspond to the inclusive event sample (see text) and 
the lower curve after the cuts x > 0.03, y < 0.3 made to suppress the pQCD contribution. (Note, 
n > corresponds to the proton beam direction.) 
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in order to get quantitative results for direct experimental considerations. 

It is of obvious interest to compare the different charm production mechanisms, which can 
easily be made through our different Monte Carlo simulation programs as discussed above. Thus, 
we use the add-on to Lepto to simulate the indirect IC scattering process, Lepto for direct 
scattering on an intrinsic charm quark and Aroma for the pQCD j*g —* cc process. In all cases 
ep events at the HERA energy were simulated with Q 2 > A GeV 2 and W 2 > 100 GeV 2 . 

In Fig.|| we display the resulting distributions in pseudo-rapidity for muons from charm 
decays in the different cases. As can be seen, the conventional pQCD process dominates strongly 
the overall rate, but can be substantially reduced by the simple cuts x > 0.03, y < 0.3 which 
hardly affects the IC contributions. The two IC processes give charm at large rapidities in the 
proton beam direction, as expected. 

One may notice that the muon detectors in the HERA experiments only covers pseudo- 
rapidities up to 3 — 4. A possibility would be to lower the proton beam energy and thereby 
decrease the strong forward boost effect. For example, lowering the proton beam energy to 
300 GeV, would essentially shift the distributions lower in rapidity by about one unit. This has 
been investigated for the direct IC scattering in ref. [pX(| , where also the background of muons 
from light meson decays (ir, K) was studied. 

Although it does not seems possible to detect these IC processes with todays set-ups in the 
HI and ZEUS experiments, it is of interest for the upgrades under discussion, e.g. within the 
presently running workshop on 'Future Physics at HERA'. 



3 Intrinsic charm in hadron interactions 

We now turn to hadron collisions for which the intrinsic charm hypothesis was first developed jlj]. 
Here also, one could consider to probe the intrinsic charm quarks in hard scattering processes 
such that pQCD is applicable. This would, however, result in a very low cross section; suppressed 
both by the IC probability and the smallness of the perturbative cross section. The charm quark 
may, however, also be put on shell and emerge in real charmed particles through softer non- 
perturbative interactions that are not strongly suppressed y|. As in previous investigations of 
this subject 0, ||], we will consider the shapes of xp distributions as derived from the IC model 
separately from the overall normalisation of the cross section. The energy dependence of the 
cross section is a non-trivial issue, which we discuss in connection with pQCD charm production 
and charm hadroproduction data. 



3.1 Basic x F distributions in IC 

A non-perturbative hadronic interaction of a proton in an IC state is thus assumed to realise 
charmed hadrons. This may occur through conventional hadronisation in a similar way as when 
the charm quark and antiquark is produced in a hard interaction, e.g. in a hadron-hadron or 
e + e~ collision. This can be well describe by a fragmentation model, like the Lund model [O], or 



a parametrised fragmentation function such as the Peterson function [18|. The effective charm 
fragmentation function is quite hard such that a simple approximation is to use a 5-function to 
let the charmed hadron get the same momentum as the charmed quark. Given the uncertainties 
of the intrinsic charm momentum distribution to be used as input for the fragmentation, this 
(5-function fragmentation is adequate here and has also been used in previous IC studies Q . 

The charm (anti) quark may also coalesce with another quark or diquark from the same 
initial proton state \uudcc) to form a hadron. This may happen when such a parton system 
forms a colour singlet, but has too small invariant mass to hadronise according to the Lund 
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model. Such a mechanism is especially important for the production of A c , but should also 
play an important role in the production of D-mesons. For an initial proton, this gives an 
asymmetry between D mesons, which can be formed by this coalescence process, and D mesons 
which cannot. Such differences in leading particle spectra are observed experimentally (see || 
and references therein) and may be an indication for these mechanisms. The c and c can also 
coalesce to form a charmonium state, mainly J/i/j or some higher resonance that mostly decays 
into the lowest lying J /tp state. 

The relative probabilities for these processes are not known, but have been discussed before. 
Following |Q, H we use the recombination probabilities 50% to form a -D-meson and 30% for a 
A c . The probability to directly form a J/ip (i.e. the cc pair is combined) is taken to be 1%. The 
absolute rates of A c and J/tp do, of course, depend strongly on these values. The D-meson rates 
and distributions are, however, not so sensitive to these values because the resulting distributions 
from fragmentation and coalescence are rather similar. The core quarks that do not coalesce 
with spectator partons, are hadronised to D-mesons with the mentioned 5-function. This is the 
case for the remaining c or c-quark in events with coalescence and both c and c in the remaining 
19% of events without coalescence. 

The momentum of the hadron formed through coalescence is taken as the sum of the corre- 
sponding ccj's, e.g. x\ c = x c + x u + Xd- The momentum distribution is then obtained by folding 
Eq. (|2|) with the proper 8 function, e.g. 5(xA c — x c — x u — Xd), and integrating out all extra degrees 
of freedom. When the core quarks are hadronised with the 5-function, the D-meson takes the 
whole intrinsic charm quark momentum as given by Eq. (||). This procedure is consistent with 
low-pt charm hadroproduction data ||. 

Based on this model we then obtain the following correlated probability distributions in 
longitudinal momentum fraction (Feynman-x) for the different combinations of charmed hadrons 
in an event; a D D pair produced through fragmentation 

dP a 180Q (l-x D -x D ) 2 x 2 D x 2 5 _ 
dx D dxp (x D + x^) 2 

a D from fragmentation and a D from coalescence 

- dP , oc 3600 (l-x D - x 5 )x 2 D (x D + x D ^ 2x 2 D h Xd \ X A ; (11) 

dx D dxQ \ x D + Xf) x d J 

a D from fragmentation and a A c from coalescence 

dP oc 36004 (^ + 3x D x Ac -x D (3x 3 + 2x A J In x d + Xk A ■ (12) 
ax A c ax £ \ l x D J 

and a J/f from coalescence 

dP 

k20(1-xj^) 2 x 3 j/iP . (13) 

ax J/tp 

These two-dimensional distributions are plotted in Fig.|(| Each one is here normalised to unit 
integral and their relative weight are the percentages discussed. Obviously, the sum of the two 
x's in each case cannot exceed unity and the distributions therefore vanish when passing the 
diagonal in the two-dimensional plots. The J/i/j, formed by the coalescence of the cc-pair, take 
a larger momentum fraction since it is formed of two quarks with relatively high xjr, as shown 
in Fig.|(| The A c , which coalesce with two valence quarks, gets a larger momentum than the D 
mesons as demonstrated in Fig.^. 
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Figure 6: Probability distributions dP j dx\dx<i and dP/dxjj^ in longitudinal momentum frac- 
tions (xi) for charmed particles from IC: (a) D D from fragmentation, Eq. (fT^j. (b) D from 
fragmentation and a D from coalescence, Eq. (11). (c) D from fragmentation and A c from 
coalescence, Eq. (li). (d) J/tp from coalescence, Eq. ftTl 



In addition to these charmed particles, higher mass states may also be produced. How- 
ever, these decay rapidly to the treated ones and can therefore be considered included in 
these parameterisations. D s cannot be produced in the coalescence process, but is allowed in 
the fragmentation process although at a suppressed rate {e.g. by a factor ~ 7 in the Lund 
model). Since its spectrum would be essentially the same as the other D-mesons, we in- 
clude it with them and use D as a generic notation for all pseudoscalar charm mesons. The 
vector mesons D* decays strongly to L>-mesons, but the decays are not charge symmetric; 
D *o D°{W0%),D* ± D°{68%)/D ± {32%). This effect is taken into account after having 
chosen a primary D* or D from spin statistics. 

The shapes of the Xir-distributions for charmed particles from intrinsic charm are thereby 
specified and found to be quite hard, as expected. In fact, they are harder than those for charm 
from perturbatively produced charm, as shown below. This gives a possibility that IC may 
contribute substantially or even dominate at large xf, depending on the absolute normalisation 
of the production cross section. 
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Figure 7: Energy dependence of charm production cross section in pp(pp) ■ The experimental data 
points are 1:^] 3:$$, 4,4':^, 5:f^, 6:%2$, 7: $2^, 8:^J 9:§F$, 10:^, 11:^ |Z 



1B:[S1, SCj and 13:M^j. The solid line (MC) is the result of our Monte Carlo calculations 
based on conventional pQCD and hadronisation, with dotted lines from variations using a naive 
application of the MRS parametrisations G and Dq of parton densities in the proton. The 
dashed curves are from our model for intrinsic charm, with the two assumptions IC1 and IC2 
about its energy dependence. For reference, the total pp cross section a tot is also shown. 



3.2 Normalisation and energy dependence of the IC cross section 

As mentioned, the main uncertainty in the intrinsic charm model is the absolute normalization of 
the cross section and its energy dependence. Since the process is basically a soft non-perturbative 
one, it may be reasonable to assume that its energy dependence is the same as for normal inelastic 
scattering Q. We therefore take as our first case 

IG1: a JC (s) = 3-W- 5 a inel (s) (14) 

shown as curve IC1 in Fig. ^ and with normalisation from ref. ||, where the magnitude of 
the cross section was estimated from data at relatively low energies (y/s = 20 — 30 GeV). 
Alternatively, one might argue that there is a stronger energy dependence related to some 
threshold behaviour for putting the charm quarks on their mass shell. We make a crude model 
for this by taking a constant fraction of the pQCD charm cross section, where such a threshold 
factor is included, to be our second case of the intrinsic charm cross section 

IC2: a W (s) = 0.1 o- pQC d(s) (15) 

as shown by curve IC2 in Fig.[?]. This is similar to the low energy (y/s = 20 — 40 GeV) treatment 
in Q. The normalisation is fixed to be the same as IC1 at the low energy where evidence is 
claimed for intrinsic charm ||. There is, however, some indication against such an increased cross 
section, as in IC2, since no evidence for J /ip from intrinsic charm was found in an experiment 
at a somewhat higher energy (800 GeV proton beam energy). 
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In Fig. we have also compiled various data on the charm production cross section in proton- 
(anti)proton collisions at various energies. A few comments on the data in Fig. |^ are here in order. 
A given experiment is only sensitive to some channels and a limited kinematical region. The 
total charm cross section is therefore obtained by a rescaling with charm decay branching ratios 
and by using assumed shapes of the xf distributions to extrapolate to unmeasured regions. In 
particular, corrections to points 1,2,6 and 7 are small while they are large for point 9 and 13. The 
bands 8,10,11 and 12 illustrate the uncertainties in these experiments due to this extrapolation. 



In band 8 the uncertainty includes a scaling for including L^-mesons (taken from [24, ^]). Data- 
band 11 is based on D + D identification, 12 on A^D and 10 on Aj!~ identification. Furthermore, 
points 3, 4 and 5 are from beam dump experiments on heavy nuclear targets without direct 
charm identification and have an additional uncertainty from the scaling with nuclear number. 
In point 4 the scaling A 0,75 has been assumed, which we have rescaled in 4' to a ^-dependence 



in order to be consistent with the other beam dump experiments and with our model [34]. Data 
points 2 and 6 come from pp interactions with explicit charm particle identification. Although 
these issues leave some uncertainty for each individual result, the combination of all data should 
give a trustworthy knowledge on the charm cross section and its energy dependence. 

These inclusive charm cross section data can be reasonably well understood by pQCD and 
does not leave much room for an IC component. This follows from our detailed Monte Carlo 
study |54j of pQCD charm production in high energy hadron collisions based on the Pythia 
program |l3| and applied to high energy cosmic ray interactions. The calculation is based on 
the conventional folding of the parton densities in the colliding hadrons and the leading order 
QCD matrix elements, for the gluon-gluon fusion process gg — > cc and the quark-antiquark 
annihilation process qq — > cc. Thus, the cross section is 

a = 1 1 1 dxidx2di /i^i'^) h( x ^Q 2 ) J (16) 

where x% are the parton longitudinal momentum fractions in the hadrons and the factorisation 
scale is taken as Q 2 = {m\ c + m 2 L - E )/2. The parton level pQCD cross section a depends on 
the Mandelstam momentum transfer t. Next-to-leading-order corrections are known and give 
a correction which can be approximately taken into account by an overall factor K = 2. The 
charm quark mass threshold s = x\X2S > Am 2 is important and is fully included in the matrix 
elements. The dominating contribution to the cross section comes from the region close to this 
threshold, since da/ds is a steeply falling distribution. It is therefore important to use QCD 
matrix elements with the charm quark mass explicitly included. The numerical value used is 
m c = 1.35 CeF/ c 2 together with A QCD = 0.25 GeV (from MRSG |§). 

The results of this pQCD calculation is also shown in Fig.]?]. At the highest energies, the 
parton densities are probed down to x ~ 10~ 5 or even below. The recent data from HERA [36, [3?]] 



show a significant increase at small x, xf(x) ~ x~ a and constrain the parton densities down to 
x ~ 10~ 4 . These data have been used in the parameterisation MRSG |35| of parton densities 
resulting in the small- x behaviour given by the power a = 0.07 for sea quarks and a = 0.30 for 
gluons. This is the most recent parameterisation, using essentially all relevant experimental data 
and can be taken as a standard choice. The effect on the total charm production cross section 
from the choice of parton density parameterisation is illustrated in Fig.|7[ The result is shown 
with MRS Dq |38| with small-x behavior x f(x) ~ const, which before the HERA data was an 
acceptable parameterisation. At high energy there is a large dependence on the choice of parton 
density functions. The difference between the G and the Dq parameterisations should however 
not be taken as a theoretical uncertainty. First of all the Do parameterisation is known to be 
significantly below the small-x HERA data and gives therefore a significant underestimate at 
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large energies. Secondly, the naive extrapolation of the G parameterisation below the measured 
region x^lO -4 at rather small Q 2 (~ m 2 ) leads to an overestimate. A flatter dependence like 
x~ e with e ~ 0.08 as x — > can be motivated ((3^] and references therein) based on a connection 
to the high energy behaviour of cross sections in the Regge framework. The implementation of 
this approach in Pythia makes a smooth transition to this dependence such that the parton 
densities are substantially lowered for x^lO -4 leading to a substantial reduction of the charm 
cross section at large energies, as given by the solid curve in Fig.[7|. 

The pQCD calculation gives a quite decent agreement with experimental charm production 
data over a wide range of energies. Nevertheless, the uncertainties in the data and the calcula- 
tions cannot exclude some smaller non-perturbative contribution. Charm production in pQCD 
is theoretically well defined and has only some limited numerical uncertainty due to parameter 
values and NLO corrections which, however, can be examined and controlled. Non-perturbative 
contributions to charm production are, however, not theoretically well defined due to the gen- 
eral problems of non-perturbative QCD. It is therefore, reasonable to take pQCD as the main 
source of charm and consider, e.g. , intrinsic charm as an additional contribution. We show 
this contribution in Fig. based on the two assumed energy dependences IC1 and IC2. Both 
cases can be consistent with the indications for IC mentioned in the Introduction, except the 
data in [^] which disfavours the stronger energy dependence of IC2. One should note in this 
context, that an even stronger energy dependence would be needed if the old ISR data were to 
be interpreted with intrinsic charm as the dominating source. 

3.3 Charm distributions at the Tevatron and LHC 

Given the differences in the pQCD and the IC mechanisms, one expects characteristic differences 
in the spectra of produced charmed hadrons at collider energies. Charm produced through the 
pQCD mechanisms should emerge with rather small longitudinal momentum or xf- This results 
from the parton fusion being largest close to threshold s = sx\X2 ~ 4m 2 .. In contrast, intrinsic 
charm is giving rise to charm particles at large fractional momenta relative to the beam particles, 
as explained before. The latter process may therefore dominate at some large xj?, with the cross- 
over point depending on the relative normalisations of the cross sections for the two processes. 

To study this, we have used our models and calculated the spectra in xp and rapidity of 
charmed hadrons at the Tevatron and LHC shown in Figs. || and |9[ Only one hemisphere 
(xf > 0, y > 0) is shown, corresponding to the proton beam direction at the Tevatron, such 
that the other hemisphere is obtained with charge conjugation symmertry. Since the pQCD 
processes are fully Monte Carlo simulated one can easily extract the xf and rapidity y = 
ln{(E + p z )/(E — p z )}- Our model for intrinsic charm is analytic and formulated in terms of 
^-dependencies. The rapidity is then calculated using y ~ In (2xFP/m±), with P the beam 
momentum and m± the transverse mass of the charmed hadron. As usual in IC models, we 
neglect the transverse momentum although it may be expected to be of the same order as the 
charm quark mass. Including p± fluctuations of this magnitude, would only cause shifts to lower 
rapidity of about ln{^/2) ~ 0.35 which is significantly smaller than the widths of the rapidity 
peaks in Fig. |9| and therefore not change the results significantly. 

In comparison to the Tevatron, the charm cross section at LHC is larger by a factor four 
for pQCD and IC2 and a factor two for IC1. The IC distributions essentially scale in xf, but 
shift to larger rapidity at LHC due to the higher beam momentum. The increase for pQCD is 
mainly at small xf, due to the charm threshold moving to smaller momentum fractions x in the 
colliding particles. 

With the stronger energy dependence of IC2, the IC cross section is significantly higher than 
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that from pQCD at large longitudinal momentum. The milder energy dependence of IC1 gives 
a charm cross section which is in general much smaller than the pQCD one and cannot really 
compete even at high xp. The only exception is for D where the previously discussed leading 
particle effects are important, resulting in cross sections of similar magnitude as pQCD for very 
large xf or rapidity. Thus, it is will be very hard to test the IC1 case, but the IC2 case could 
be observable provided that the forward coverage of the detectors is extended far enough. This 
could be considered in connection with dedicated heavy flavour experiments covering the forward 
region in particular. 



4 Summary and conclusions 

The hypothesis of intrinsic charm quark-antiquark pairs in the proton wave function is not 
ruled out by experiment. On the contrary, there are some evidence in favour of it, but no safe 
conclusion can be made at present. It is therefore important to consider various ways to gain 
additional information that could help clarify the situation. Based on previous work, we have 
constructed explicit models for intrinsic charm and how it may be examined in deep inelastic 
lepton-nucleon scattering and hadron collisions. Our models are partly implemented in terms of 
Monte Carlo programs, which allow detailed information to be extracted since complete events 
are simulated. 
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Figure 8: Cross section versus longitudinal momentum fraction xf for charm particles produced 
from pQCD (full curves) and from the intrinsic charm model with energy dependence IC1 (dashed 
curves) and IC2 (dotted curves). In each case, the upper curves correspond to LHC (pp at 
yfs = lATeV) and the lower ones to the Tevatron (pp at l.STeV). Only one hemisphere 
xf > (proton beam direction) is given. 
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Figure 9: Cross section versus rapidity for charm particles produced from pQCD (full curves) 
and from the intrinsic charm model with energy dependence IC1 ( dashed curves ) and IC2 ( dotted 
curves). In each case, the upper curves correspond to LHC (pp at y/s = \ATeV) and the lower 
ones to the Tevatron (pp at 1.8TeV). Only one hemisphere y > (proton beam direction) is 
given. 



In DIS such intrinsic charm quarks can be directly probed and we find that it may dominate 
the inclusive charm F| structure function at large x. Muons from charm decay can be used to 
signal events with charm and we compare results of the IC mechanism with the conventional per- 
turbative QCD boson-gluon fusion process, as well as background muons from n and K decays. 
We devise cuts in order to enhance IC relative to the backgrounds. Signal and backgrounds can 
then be brought to about the same level. We point out that data samples already collected by 
NMC would be suitable for this purpose and could be sensitive down to a level of about 0.2% 
probability of intrinsic charm in the proton. 

Direct scattering on intrinsic charm quarks at HERA has been investigated before |p.Q|. 
Here, we investigate the indirect process with scattering on a light valence quark such that 
the cc fluctuation in the proton remnant is realised. The rapidities of the produced charmed 
particles have been calculated and found to be very forward, as expected. Present detectors 
do not have enough forward coverage to detect these processes, but one may consider them in 
connection with possible upgrades for future HERA running. 

For hadronic interactions, the intrinsic charm model gives definite predictions for charmed 
particle xp spectra. The absolute normalisation and its energy dependence is, however, not 
clear. We have investigated this in comparison with conventional pQCD productions mechanism 
for charm and measured charm production cross sections. This constrains the allowed energy 
variation of the IC cross section. Using two simple models for this energy dependence, we 
calculate the xp and rapidity distributions for charmed particles at the Tevatron and LHC. We 
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find that it is only if the IC cross section has a significant rise with energy, that it can compete 
with normal pQCD production of charm. In any case, the IC contribution is at very large 
forward momenta, such that its detection would require coverage at very forward rapidities. 

In the context of very high energy hadronic collisions one should also consider cosmic ray 
interactions in the atmosphere, which we have investigated in detail [34]. Here also, the intrinstic 
charm mechanism would contribute significantly only if it has a strong energy dependence. 

Finally, one should remember the possible extension from intrinsic charm to intrinsic bottom 
quarks. Although we have not presented numerical results, our methods can easily be applied for 
intrinsic bottom processes. In comparison to charm the general expectations are, as discussed 
in the Introduction, that the ^-distributions will be only slightly harder but the overall rates 
lower by about a factor ten. 
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